Identification of an RVB liquid phase in a quantum dimer model with competing 

kinetic terms 
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Starting from the mean-field solution of a spin-orbital model of LiNi02, we derive an effective 
quantum dimer model (QDM) that lives on the triangular lattice and contains kinetic terms acting 
on 4-site plaquettes and 6-site loops. Using exact diagonalizations and Green's function Monte Carlo 
simulations, we show that the competition between these kinetic terms leads to a resonating valence 
bond (RVB) state for a finite range of parameters. We also show that this RVB phase is connected 
to the RVB phase identified in the Rokhsar-Kivelson model on the same lattice in the context of a 
generalized model that contains both the 6-site loops and a nearest-neighbor dimer repulsion. These 
results suggest that the occurrence of an RVB phase is a generic feature of QDM with competing 
interactions. 

PACS numbers: 75.10.Jm, 05.50.+q, 05.30.-d 



I. INTRODUCTION 

After their first derivation by Rokhsar and Kivelson in 
1988 in the context of cupratesji the hard-core quantum 
dimer models (QDM) have attracted significant atten- 
tion. The phase diagrams of the QDM on the square 
and triangular lattices have been investigated in great 
details and, following the pioneering work of Mocss- 
ner and Sondhi on the triangular lattice^ the very exis- 
tence of a stable resonating valence bond (RVB) phase 
has been unambiguously demonstrated^ The presence of 
a liquid phase with deconfined vison excitations^ has also 
been established for a toy model living on the Kagome 
latticed 

However, the relationship between QDMs and Mott 
insulators, the physical systems for which they were pro- 
posed in the first place, is not straightforward. It is well 
established by now that the ground state of the S=l/2 
Heisenberg model on the square and triangular lattices 
exhibits long-range magnetic order of Neel and 120 de- 
gree type respectively, and this type of order cannot 
be reached within the variational basis of Rokhsar and 
Kivelson, which consists of short-range singlet dimers. 
For a QDM to be a good effective model, one should thus 
identify models for which the subspacc of short-range 
dimer coverings on a certain lattice is a good variational 
basis. 

The first example of such a case was provided by the 
S=l/2 Heisenberg model on the trimerized Kagome lat- 
ticed Indeed, an effective spin-chirality model living on 
a triangular lattice can be derived, and, at the level 
of a mean-field decoupling between spin and chirality, 
the ground state manifold consists of all dimer coverings 
on the triangular lattice. Going beyond mean-field is 
thus expected to lead to a relevant effective QDM. Us- 
ing Rokhsar and Kivelson's prescription, which consists 
in truncating the Hamiltonian and inverting the overlap 
matrix within the basis of dimer coverings, Zhitomirsky 



has derived such an effective Hamiltonian^ and shown 
that the main competition is between kinetic terms in- 
volving loops of length 4 and 6 respectively, and not a 
competition between a kinetic and a potential term, as 
for the QDM derived by Rokhsar and Kivelson. The next 
logical step would be to study the properties of this effec- 
tive QDM. This is far from easy however. We know from 
the experience with the standard QDM model on the 
triangular lattice that the clusters reachable with exact 
diagonalizations are much too small to allow any signifi- 
cant conclusion regarding the presence of an RVB phase, 
and since there is no convention leading only to negative 
off-diagonal matrix elements, it is impossible to perform 
quantum Monte Carlo simulations. 

Recently, two of us came across another model, for 
which the low energy sector consists of almost degenerate 
singlet coverings on the triangular lattice. This model is a 
Kugel-KhomskiilS, model that was derived in the context 
of LiNiC>2, and the mean- field equations that describe 
the decoupling of the spin and orbital degrees of free- 
dom possess an infinite number of locally stable solutions. 
These solutions are almost degenerate and correspond to 
spin singlet (and orbital triplet) dimers on the triangular 
latticed Following Rokhsar and Kivelson's prescription, 
an effective QDM can also be derived (see Appendix 0. 
As for the S=l/2 Heisenberg model on the trimerized 
Kagome lattice, it consists of a competition between ki- 
netic terms, with two important differences however. The 
main term of length 6 lives on loops that have a shape 
of large triangles, a term absent in the other case. But 
more importantly, the off-diagonal matrix elements are 
all negative. 

Since the competition between kinetic processes was 
never investigated before, we have decided to concentrate 
on the minimal model obtained by keeping only the dom- 
inant term of length 6 for clarity. We have checked that 
the properties of the complete effective model are similar. 
This minimal model is described by the Hamiltonian: 
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where the sums run over the 4-site and 6-site loops with 
all possible orientations. Although the repulsion is a 
higher order process, we have included a repulsion term 
in the Hamiltonian, and we will treat its amplitude V 
as a free parameter to be able to make contact with the 
Rokhsar-Kivelson model on the triangular lattice. The 
hopping amplitudes t and t' are negative, and although 
the ratio if Jt is in principle fixed by their expression in 
the perturbative expansion, we will also treat it as a free 
parameter. 

Our central goal in this paper is to determine the na- 
ture of the ground state as a function of t' /t. With re- 
spect to what we already know about QDMs, the main 
question is whether a competition between kinetic terms 
can also lead to a liquid phase. As we shall see, the an- 
swer to that question is positive, a liquid phase being 
present in a finite region of the phase diagram in the 
t'— V plane. 

The paper is organized as follows. In Section [HI we 
briefly review the basic preliminaries used in the rest of 
the paper. The results obtained with exact diagonaliza- 
tions are presented in Section II I II those obtained with 
quantum Monte Carlo in Section llVl and the conclusions 
in Section [V] The perturbation calculation which has 
motivated the investigation of this QDM is finally pre- 
sented in Appendix lAl 



II. THE METHOD 

In this section, we present a brief introduction to the 
numerical methods, to the clusters used in the analysis, 
and to the physical concepts underlying the determina- 
tion of the phase diagram. More details can be found in 
Ref . U 

Let us first discuss the shape of the finite-size clusters. 
In general, a finite cluster is defined by two vectors Ti 
and T2 and, in order to have the symmetries for rotations 
by 27r/3, they have to satisfy^ 

Ti = Zui + TOU2 

T 2 = — mui + (I + m)u 2 , 

where I and m are integers and Ui = (1,0) and u 2 = 
(l/2,\/3/2) are the unitary vectors defining the trian- 
gular lattice. The number of sites in the cluster is 
N = I 2 + m 2 + lm. In order to have also the axial sym- 
metry, and therefore all the symmetries of the infinite 
lattice, we must take either lm = or I = m. The first 
possibility corresponds to type-A clusters (with the no- 
tation of Ref. 0), with N = I 2 sites; the second one gives 
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FIG. 1: Example of type-A (left) and type-B (right) clusters 
with 16 and 12 sites, respectively. 



rise to type-B clusters, with N = 3 x I x I (for examples 
of both cases, see Fig.^l. Since for t'/t = and V/t = 
the ground state belongs to a crystalline phase with a 
12-site unit cell^^ we will restrict in this work to clus- 
ters with a number of sites multiple of 12, in order not to 
frustrate this order. To limit the finite size effects related 
to the geometry of the clusters, we will concentrate on 
type-B clusters, which are always compatible with this 
order. Note that the 6-dimer loop kinetic term does not 
introduce further restrictions since all it requires is to be 
able to accommodate 6-site unit cells. 

A very important concept in the QDM is the existence 
of topological sectors. Indeed, in the triangular lattice, 
the Hilbert space is split into four disconnected topologi- 
cal sectors on a torus defined by the parity of the number 
of dimcrs cutting two lines that go around the two axes of 
the torus and denoted by (p, q) with p, q = (respectively 
1) if this number is even (respectively odd). One can con- 
vince oneself, by a direct inspection of the effect of the 
4-site and 6-site terms, that these numbers are conserved 
quantities under the action of the Hamiltonian of Eq. J^. 
More generally, this is a consequence of the fact that the 
topological sectors are not coupled by any local pertur- 
bation. These topological sectors are extremely useful to 
distinguish between valence bond solids and spin liquids. 
Indeed, valence bond solids are only consistent with some 
topological sectors, whereas RVB spin-liquid phases are 
characterized by topological degeneracy. Therefore, the 
main goal will be to investigate whether, in the thermo- 
dynamic limit, the topological sectors are degenerate or 
not. In that respect, it is useful to remember that the 
(0,1) and (1,0) sectors are always degenerate with ei- 
ther (0,0) or (1, 1) (depending on the cluster geometry) 
since they contain the same configurations rotated by an 
angle ir/3~ One can thus, without any loss of general- 



ity, restrict oneself to the analysis of the (0, 0) and (1,1) 
sectors. Therefore, we define the absolute value of the 
topological gap as: 

AE = \E 00 — Eix\, (2) 

where Eqq and En are the total ground-state energies 
for the topological sectors with p = q = and p = q = 1 , 
respectively. This gap is expected to scale to zero with 
the cluster size in the RVB phase. 

Finally, in order to detect a possible dimer order, we 
also consider the static dimcr-dimcr correlations 

D^ir-r') = (D l (r)D3(r')), (3) 

where D l (r) is the dimer operator defined as follows: It is 
a diagonal operator in the configuration space that equals 
1 if there is a dimer from the site r to the site r + ai, with 
oi = (1,0), a 2 = (1/2,^3/2), or a 3 = (-1/2,^/2) and 
vanishes otherwise. 

The method used is the same one as that used by Ralko 
and collaborators^ to determine the phase diagram of the 
Rokhsar-Kivelson QDM on the triangular lattice and our 
investigations are Lanczos diagonalizations and Green's 
function Monte Carlo (GFMC) simulations. In particu- 
lar, the GFMC is a zero-temperature stochastic technique 
based on the power method: Starting from a given wave 
function and by applying powers of the Hamiltonian, the 
ground state is statistically sampled to extract its energy 
and equal-time correlation functions. In principle, as in 
other Monte Carlo algorithms, in order to reduce the 
statistical fluctuations, it is very useful to consider the 
importance sampling, through the definition of a suit- 
able guiding function. Unfortunately, when dealing with 
dimer models, it is very hard to implement an accurate 
and, at the same time, efficient guiding function for the 
crystalline phases. This problem is particularly relevant 
when the 6-site term becomes dominant. In this case, 
our simulations suffer from wild statistical fluctuations, 
deteriorating the convergence of the GFMC. As a con- 
sequence, we are not able to reach the largest available 
size, 432-site cluster, for all parameters t'/t and V/t. By 
contrast, given the simple form of the spin-liquid ground 
state (that reduces to a superposition of all the config- 
urations with the same weight at the Rokshar-Kivelson 
point), in the disordered region we can use the guiding 
function with all equal weights and obtain very small fluc- 
tuations (no fluctuation at the Rokshar-Kivelson point) 
and, therefore, excellent results with zero computational 
effort. Combining these facts, the loss of the GFMC con- 
vergence can be interpreted as a signal for the appearance 
of a crystalline phase. Of course, this is not a quantita- 
tive criterion, but, as it will be shown in the following, 
it gives reasonable insight into the emergence of a dimer 
order. 

Finally it should be mentioned that, since the different 
topological sectors are completely decoupled (each dimer 
configuration belonging to one and only one of them) 
and cannot be connected by the terms contained in the 
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FIG. 2: (Color online) Difference between _Eoo and -En, the 
total ground-state energies of the topological sectors with p = 
q — and p = q = 1, respectively. The results are found by 
exact diagonalizations for clusters with 12 and 36 sites. 




FIG. 3: (Color online) Dimer-dimer correlations on the 36- 
site cluster. The dimer of reference is the thickest one in the 
up-right corner. The thicker the line, the farther the value 
of the correlation from the uniform distribution equal value 
1/36. Solid lines are used when the correlation is higher than 
1/36, and dashed lines when it is lower. 



Hamiltonian, within the GFMC it is possible to work in 
a given topological sector, making it possible to extract 
the ground-state properties of each of them. 



III. EXACT DIAGONALIZATIONS 

To get a first idea of the properties of the model, we 
start with the results we have obtained with exact di- 
agonalizations of finite clusters for the model of Eq. Q 
with V/t = 0. Let us first begin with the ground-state 
energy for the 12- and 36-site clusters for both the topo- 
logical sectors (0,0) and (1,1) (sec Fig. Note that 
the 12-site cluster is of type B, whereas the 36-site one is 
of type A. We have that, for both sizes, a level crossing 
occurs for t'/t ~ 2. Below that value, the ground state is 
in (1, 1) topological sector, in agreement with the earlier 
results of Ref. for V/t = 0. The main difference be- 
tween the 12-site and the 36-site clusters is that, for the 
36-site cluster, the topological ground-state energies stay 
very close in a large parameter range for t'/t ~ 2. This 
could suggest that, upon increasing the size, this level 
crossing might evolve into a phase where these energies 
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are rigorously degenerate, giving rise to a liquid phase 
without any crystalline order. 

In order to give an idea of the various phases, we re- 
port in Fig. the dimcr-dimer correlations of Eq. Q 
for the 36-site cluster below, at, and above the level 
crossing, which takes place at t'/t = 2.2 for this size. 
For small t'/t, the correlations show a pattern similar to 
that of the intermediate \/Vl x ^/l2 phase of the stan- 
dard QDM model, already shown in Ref. 0- It should 
be stressed that, since in these calculations the Hamil- 
tonian has the translational symmetry, the 12-site unit 
cell is not directly visible from Fig. and in order to 
have a clearer evidence one should break the symmetry 
by hand. Nonetheless, as it has been shown in Ref. Oj 
these results are in perfect agreement with the existence 
of a \JV1 x \JY1 phase with a crystalline ground state 
in the thermodynamic limit. When t'/t is large, another 
pattern arises, which has never been observed in the stan- 
dard QDM, and which presents a kind of 6-site triangle 
ordering. In this case, the dominant kinetic term involv- 
ing 6 sites [see the second term of the Hamiltonian 
induces a dimer pattern with the same symmetry, pos- 
sibly inducing a ground state with a 6-site unit cell in 
the thermodynamic limit. Also in this case the trans- 
lational symmetry of the ground state partially masks 
the existence of a regular dimer pattern. Unfortunately, 
we will not be able to confirm this prediction since, as 
stated before, the GFMC algorithm has serious problems 
of convergence inside this phase and for large clusters. 

In any case, for intermediate values of t'/t, the correla- 
tions decay very rapidly with the distance and are close to 
those obtained in the liquid phase of the standard QDM 
with V/t < 1 and t'/t = 0£ This fact gives another evi- 
dence of the possible existence of an RVB phase between 
two ordered phases in the model with competing kinetic 
terms and without the dimer repulsion. Of course by 
considering the exact diagonalization results only it is im- 
possible to give definite statements on the stabilization of 
this liquid phase and, therefore, in the following section, 
we will consider a more systematic study of the topolog- 
ical gap, in order to unveil the existence of a wide dis- 
ordered region that develops from the Rokshar-Kivelson 
point V/t = 1 of the standard QDM and survives up to 
V/t = and finite t'/t. 



IV. GREEN'S FUNCTION MONTE CARLO 

In this section, we use the GFMC method to extend 
the results of the previous section to larger clusters, with 
up to 432 sites, and to map out the phase diagram in the 
t'—V plane. 



A. The case of V/t = 

Let us first describe the results we have obtained for 
V/t = and consider the behavior of the topological gap 
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FIG. 4: (Color online) Topological gaps for V/t — as a 
function of 1/v^V, where N is the number of sites and for 
different values of t'/t. Upper panel: Small values of t'/t, 
where the gap opens for large clusters. Lower panel: Larger 
values of t'/t. For t'/t > 1.6, the finite-size gap closes upon 
increasing the cluster size, signaling a liquid phase. 



given in Eq. J3J). In Fig. 01 for clarity, we divided the re- 
sults in two sets for < t'/t < 0.6 and t'/t > 1. The first 
remarkable feature is that, for most parameters, the gap 
decreases between 12 and 48 sites, regardless of its be- 
havior for larger sizes. Therefore, the possibility to study 
much larger clusters is crucial for this analysis. Indeed, 
there is a clear change of behavior for t'/t ~ 1.6, a value 
beyond which our extrapolations give solid evidences in 
favor of a vanishing gap in the thermodynamic limit. On 
the other hand, for smaller t'/t ratios, we have a clear 
evidence that AE increases for large sizes. Therefore, we 
come to the important conclusion that, the crystalline 
Vl2 x a/12 phase is destroyed by increasing the ampli- 
tude of the 6-site term and the system is eventually driven 
into a liquid RVB phase. 

To have a further confirmation of the existence of this 
disordered phase, we have calculated the dimer-dimcr 
correlations. The results obtained on the 108-site cluster 
for t'/t = 2 are reported in Fig. |SJ where the correla- 
tion functions for parallel dimers along the same row are 
plotted as a function of the distance. Given the small 
number of clusters available, a precise size scaling of the 
order parameter is not possible and also a meaningful 
estimation of the correlation length is very hard. Nev- 
ertheless, the behavior is definitely consistent with an 
exponential decay and the uncorrelated value of 1/36 is 
approached very rapidly, as expected in a liquid phase 
without any crystalline order. 

Unfortunately, the larger t'/t region is numerically far 
more difficult to access. Indeed, as stated above, al- 
though the GFMC is in principle numerically exact, we 




FIG. 5: (Color online) Dimer-dimer correlation function for a 
108-site cluster along a horizontal line as a function of the dis- 
tance for t'/t = 2 and V/t = 0. The dashed line corresponds 
to 1/36, the value in the absence of correlations. 



have not been able to find an efficient guiding function 
to perform the importance sampling and wild statisti- 
cal fluctuations prevent us to reach a safe convergence 
for large clusters. In practice, we have access to clusters 
up to 108 sites, that are still too small to predict the 
thermodynamic behavior. For instance, for V/t = 0, the 
convergence stops before one can observe any increase of 
the topological gap, and the criterion used for the phase 
transition on the other side of the RVB phase cannot 
be used any more. However, the lack of convergence is 
a clear sign that one enters a new (crystalline) phase. 
So, if the change of behavior of the topological gap with 
the size cannot be observed, we take as a definition of 
the boundary for the phase transition the parameters for 
which the convergence is not good any more. This is ex- 
pected to be semi-quantitative, and indeed, as we shall 
see in the next section when studying the full diagram, 
the points obtained with this criterion agree reasonably 
well with those obtained with the reopening of the topo- 
logical gap. 

In summary, although the region where the 6-site term 
dominates over the usual 4-site dimcr flip is not accessible 
by using GFMC, based on our numerical results for small 
t'/t, we can safely argue that the crystalline \/l2 x y/l2 
phase is destabilized by increasing the 6-site kinetic term, 
leading to a true disordered ground state with topological 
degeneracy. 



B. Phase diagram in the t'—V plane 

In this section we prove that the RVB phase found in 
the previous paragraph for V/t = is connected to the 
one obtained for the standard QDM, i.e., close to the 
Rokshar-Kivelson point and t 1 ft = 0. In order to do 
that, we have investigated a generalization of the pre- 
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FIG. 6: (Color online) Topological gap for t'/t = 1 and for 
various V/t as a function of 1/viV, where N is the number 
of sites. Small and large values of V/t have been shown sep- 
arately in the upper and lower panels for clarity. 



vious model that also includes a repulsion V between 
dimers facing each other, sec Eq. Q). For t'/t = 0, this 
model reduces to the standard QDM, which has an RVB 
phase for 0.75 < V/t < To map out the complete 
phase diagram of this model, we have done the same anal- 
ysis of the previous paragraph for different values of V/t 
between and 1. As an example, we show in Fig. E3 the 
finite-size scaling of the topological gap for t'/t = 1 and 
several values of V/t. The global behavior is the same 
as before and we have clear evidence that the topological 
gap present at V/t = persists up to V/t ~ 0.25, and 
that it opens again for V/t ~ 0.8. In that case, as in 
many other cases, it turned out to be possible to actu- 
ally observe the opening of the gap upon leaving again 
the RVB phase before the convergence problems were too 
strong. 

The resulting phase diagram is depicted in Fig. \7\ 
where different symbols have been used depending on 
whether the boundary was determined from the increase 
of the correlation length with the size or from the loss of 
convergence: Solid circles when the closing of the topo- 
logical gap was observable and solid diamonds when the 
convergence of the GFMC was lost for large systems. In- 
terestingly enough, these different criteria build a rela- 
tively smooth line, a good indication that it can be in- 
terpreted as a phase boundary. Note that we did not 
perform simulations for V/t > 1 where, for a vanishing 
t'/t, a crystalline phase with staggered dimer order is 
stabilized. So we cannot exclude that the boundary ex- 
tends beyond V/t > 1 for small t'/t. Remarkably, the 
RVB phase we found for V/t = is connected to the 
RVB phase reported before for the standard QDM, i.e., 
t'/t = 0, the total RVB phase building up a large stripe 
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FIG. 7: (Color online) Phase diagram in the t' — V plane. A 
wide disordered region extends all the way from the standard 
QDM (t'/t = axis) to the purely kinetic QDM (V/t = 
axis). The description of the symbols is given in the text. 

that encompasses a significant portion of the phase dia- 
gram. We have also calculated static correlation function 
for several values of the parameter, but they merely con- 
firm the identification of the phases and are not reported 
for brevity. 

V. CONCLUSIONS 

Coming back to our long-term motivation, namely to 
find an RVB phase in a realistic model of Mott insula- 
tors, this paper contains significant results of two sorts. 
First of all, we have shown (see Appendix |A"|) that, start- 
ing from the quasi-degenerate mean-field ground state 
of a Kugcl-Khomskii spin-orbital model, one can con- 
struct a QDM with two remarkable properties: It de- 
scribes a competition between two kinetic terms of com- 
parable magnitude, and all off-diagonal matrix elements 
in the dimer basis are negative. This has allowed us to 
implement the GFMC and to investigate the results of 
the competition between these terms. It turns out that 
the competition between these kinetic terms leads to the 
disappearance of the \/l2 x VT2 crystalline order when 
t'/t ~ 1.6. This transition is similar to the transition into 
the RVB phase that happens in the standard QDM upon 
approaching the Rokhsar-Kivclson point. Indeed, the two 
phases can be connected into a single RVB phase in the 
context of a generalized QDM. As far as the numerical in- 
vestigation of the model is concerned, the main open issue 
is to pin down the nature of the phase that occurs when 
the 6-dimer kinetic term dominates. Unfortunately, the 
GFMC suffers from severe statistical fluctuations when- 



ever the guiding function is not accurate, i.e., for clusters 
larger that 108 sites and large t'/t. Therefore, we cannot 
make any definite statements on the phase where the 6- 
site term dominates. Another interesting question is of 
course the nature of the quantum phase transitions be- 
tween these phases (continuous or first order). We are 
currently working on that rather subtle issue in the con- 
text of the standard QDM. 

The general features of the RVB phase are consistent 
with the phenomenology of LiNi02, which exhibits nei- 
ther orbital nor magnetic long-range order. A number of 
points deserve further investigation however. The pre- 
cise form of the QDM does not seem to be an issue: 
The actual model that can be derived along the Rokhsar- 
Kivelson lines has more terms (see Appendix^), but pre- 
liminary results show that the RVB phase is present in 
that model as well. The fact that the RVB phase does not 
contain the point t'/t = 1.34 derived in the Appendix lAl 
is not really an issue either. First of all, this ratio was 
determined for vanishing Hund's rule coupling and one 
vanishing hopping integral, and its precise value in that 
case should at best be taken as an indication of its order 
of magnitude in the actual system. Besides, the other 6- 
site terms pull the RVB region down to smaller values of 
the relative ratio of the 6-dimer term to the 4-dimer term. 
What would deserve more attention is the validity of the 
expansion that leads to the effective QDM. The small 
parameter of the expansion is not that small (a = l/v2, 
see Appendix |A"|> . and it would be very useful to bet- 
ter understand to which extent such an expansion can 
be controlled. Nevertheless, the present results strongly 
suggest that the presence of an RVB liquid phase be- 
tween competing ordered phases is a generic feature of 
QDM, and that to identify such a phase in a realistic 
Mott insulator via an effective QDM is very promising. 

We acknowledge useful discussions with P. Fazekas, M. 
Ferrero, and K. Pcnc. This work was supported by the 
Swiss National Fund and by MaNEP. F.B. is supported 
by CNR-INFM and MIUR (COFIN 2005). 



APPENDIX A: DERIVATION OF THE MODEL 

LiNi02 is a layered compound in which the Ni 3+ ions 
are in a low spin S=l/2 state with a two-fold orbital 
degeneracy. A fairly general description of this system is 
given by a Kugcl-Khomskii Hamiltonian defined in terms 
of two hopping integrals t^ and t' h , the on-site Coulomb 
repulsion U and the Hund's coupling J which, on a given 
bond, takes the formii 
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with the usual definitions for the projectors on the singlet 
and triplet states of a pair of spins: 
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and Vfr 1 



(A2) 



The vector n?- depends on the type of bond. With the 
convention of Fig. [H] they are given by: 



(0,0,1) 

1 l 

2 ' 

—-) 

\ 2 ' u ' 2> - 



n 12 

nf 3 = (^,0,-i) 



(A3) 



The operators are pseudo-spin operators acting on the 
orbitals. 

On a given bond (see Fig.|HI), a dimer is defined by the 
following wave function: 



(A4) 



where \4>ij) an( i Ifilj) are respectively the spin and orbital 
components. The spin component is the usual singlet 
given by: 



l^> = a(IT4j>-UiTi», 



(A5) 



regardless of the orientation of the bond. We have de- 
noted by a the normalization coefficient, whose explicit 
value is of course given by l/\/2, to be able to keep track 
of the order in a of various overlaps. The orbital part 
depends on the bond and is given by: 



\ai)\a 2 ) 



73 



:a 3 



I IT X I 1 , x I 1 

1023) = I - + — b 2 )\ - -a 3 



V3, 



^3> 

2 7 



(A6) 
(A7) 

(A8) 



with the convention of Fig. [8] and with \a) = \d^, z i- r 2) 
and | b) — \d x 2_ y 2). We also use the convention that the 
wave function has a plus sign if the two sites are in the 
order defined by the arrows of Fig. |H1 and a minus sign 
otherwise. 

With these definitions, the matrix element of the 
Hamiltonian on each bond is given by: 



e 2a 1 



(A9) 




FIG. 8: Sign convention for a triangle 




FIG. 9: A 4-dimer loop 
where e is defined by: 



< 



(A10) 

U 2 - J 2 

According to the mean-field analysis of Ref. 0, the 
wave functions obtained as tensor products of these wave 
functions on all dimer coverings on the triangular lattice 
should constitute a good variational basis if t' h <§; th and 






FIG. 10: The 3 types of 6-dimer loops 
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TABLE I: Overlap and matrix element for two configurations 
differing by a loop-4, a loop-6 of type (1), (2) and (3). 





loop-4 


loop-6 (1) 


loop-6 (2) 


loop-6 (3) 


{m\H\n) 


66ea 4 
16^ 

~ -0.0644e 


103ea B 
16 2 

~ -0.0503e 


24ea° 

~ -0.0012e 


177ea H 

~ -0.0864e 



J < [f. In the following, wc consider for simplicity the 
case t' h = and J = 0. 

Following Rokhsar and Kivelson^I the idea is now to 
perform a unitary transformation to derive the effec- 
tive QDM. If one defines the overlap matrix by S mn = 
(^ml^n), where |^ m ) and \^ n ) are dimer coverings, the 
states defined by 

H = EM) ( a11 ) 

— \ / m,n 

n 

constitute an orthonormal basis and the matrix elements 
of the Hamiltonian in this basis are given by: 

H c ± ~ (m\n\n) = J2 (S-*) , (* fc |H|*j) M), . 

kl 

(A12) 

The inverse of the square root of the overlap matrix can- 
not be calculated exactly, but this can be done approx- 
imately in the context of an expansion in powers of a. 
Indeed, the overlap matrix can be expanded as: 

S = I + 2Aa 4 + 2Ba 6 + 0[a 8 } , (A13) 



which leads to: 

{Sy? = I - Aa 4 - Ba 6 + 0[a 8 }. (A14) 

In these expressions, the matrices A and B only have non 
vanishing matrix elements, equal to 1, between configu- 
rations that are the same except on a 4-dimer loop (see 
Fig. |5J) or on one of the three types of 6-dimer loop (see 
Fig. irU|) . respectively. 

Similarly, the Hamiltonian matrix TL defined by Tl mn = 
m \H\^> n ) ~ eNdS mn , where Nd is the number of dimers, 
has an expansion in powers of a that reads: 

U = Ca 4 + Da 6 + 0[a 8 }, (A15) 

where the matrices C and D have non vanishing matrix 
elements under the same conditions as matrices A and 
B. 

Since all these expansions start with a 4 , it is clear that 
the first contribution to the diagonal part of H cS will be 
of order a 8 . So, to order a 6 , the effective Hamiltonian 
will only have off-diagonal matrix elements. Moreover, 
to this order, these matrix elements are simply given by 
(m\ii\n} = (* m |W|*„). They are tabulated in Table □ 
for configurations that are the same except on a 4-dimer 
loop (Fig.|5Jl or on one of the three types of 6-dimer loop 
(Fig.llUp. All these matrix elements are negative. Among 
the 6-dimcr loop terms, the matrix clement of type (3) 
is the largest, and its ratio to the 4-dimer loop term is 
1.34. 
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